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1 . BACKGROUND 


Serious exploration of the upper atmosphere began in the 1950's as the 
possibility of artificial earth satellites appeared on the horizon. In 1960, 
it was shown by Hines ^ that much of the observational data on the upper 
atmospheric variations could be interpreted in terms of gravity waves. Al- 
though the original theory considered only an isothermal atmosphere, with no 
viscosity or heat conduction it did explain several features which are con- 
sistent with observation. Among these features were the large ratio between 
horizontal and vertical wavelengths, the increasing amplitude with height, and 
the downward, phase propagation associated with upward energy propagation. Later 
advances in the theoretical considerations analyzed the effects of realistic 

/n \ /Q / \ 

temperature profiles^ , viscous damping and thermal conductivity^ * , wind 

shear^, critical layers^, and ohmic losses^’^. 

Although there is no unambiguous resolution of gravity waves in the upper 
atmosphere (i.e. simultaneous observation of amplitudes, phases, and frequency 
sufficient to verify propagation according to the theoretical dispersion 
equation) there is strong circumstantial evidence for their occurrence. Table 
1 summarizes some of this observational evidence. This table was adapted from 
Bauer et al v . Tables 2 and 3, adapted from Bauer's report and from Tchen v } 
summarize the mechanisms for the generation and the dissipation of gravity waves. 
In addition to the dissipation effects of Table 3, gravity waves can also be 

prevented from propagating to adjacent layers of the atmosphere by the effects 

- , _ (41) 

of ducting 

It should be emphasized, however, that despite the large body of observa- 
tions summarized in the above tables, that in no cases have sufficient parameters 


TABLE 1. EVIDENCE FOR THE OCCURRENCE OF GRAVITY 
WAVES IN THE ATMOSPHERE 


Observation 


Altitude 

Range 


References 


1. Wave structure on clouds 

2. Pressure variations observed on 
microbarographs 

3. Mountain waves (well known to 
sailplane pilots) 

4. Wave structure on wind profiles 
from ROBIN balloons and smoke 
trails 

5. Wave structure on pitot tube 
and grenade temperature pro- 
files 

6. Wave patterns on noctilucent 
clouds 

7. Irregular meteor wind varia- 
tions 

8. Chemical releases show winds 
with a wave structure 

9. Wave structure in rocket probe 
measured temperature profiles 

10. Occasional sinusoidal variation 
in ion temperature observed on 
satellite-borne instruments 

11. Thomson scatter radar obser- 
vations of waves 


0 - 10 km 


low altitudes 


0 - 10 km 


30 - 70 km 


30 - 120 km 


80 - 90 km 


80 - 100 km 


80 - 160 km 


Cole and Kantor 
Mahoney 


(9) 


Theon et al 


(ID 


Witt (12) , Hines (13 ' ) 


Haurwitz and Fogle 


(14) 


Greenhow and Neufeld 
Justus and Roper 


(15) 


Spizzichino and Revah 


Kochanski^^ 


(17) 


100 - 170 km Knudson and Sharp 


(19) 


200 - 250 km Harris, et al. 


( 20 ) 


F - region 
(200 - 300 km) Waldteufel 


Vasseur and 

( 21 ) 
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TABLE 1 . Continued 


Observation 


Altitude 

Range References 


12 . Various ionospheric irregu- 
larities, such as Sudden 
Ionospheric Disturbances 
(SIDs) 

13. Certain ionospheric disturb- 
ances correlate with severe 
local storms in the troposphere 


200 - 300 km 


F2 - region 


Georges 


(22,23) 


(24) 

Baker and Davies 
Georges 

(25) 

Davies and Jones 
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TABLE 2. POSSIBLE MECHANISMS FOR THE GENERATION 
OF GRAVITY WAVES IN THE ATMOSPHERE 

( 28 ) 

1. Generation of propagating wave by surface frontal systems (Hines' - 

fl3 28} 

2. Generation of waves induced by jet streams (Hines v 9 J* 

3. Winds blowing transversely over mountain ranges, or equivalent effects on 

the air motion produced by warm cities, etc., either directly or by vortices 

(29) 

produced by air flows over non-uniform surfaces (cf. Eckart , p. 120). 

4. Volcanic eruptions, thunderstorms, hurricanes, typhoons, tornadoes, and 
similar dramatic meteorological disturbances. 

5. Temperature inversions in the lower portion of the atmosphere produce a 

fluctuating wind, which correlates with gravity waves observed on clouds 

(29) 

at the lower levels (cf. Eckart , p. 121). 

6. % Solar heating and photodissociation with the stratosphere/mesosphere, 

leading to local expansion. 

7. Photochemical "pumping 11 by either oxygen dissociation (Leovy) or ozone 
dissociation (Shere and Eberstein ) . 

(3D (32) (33) 

8. Auroral currents (Blumen and Hendl v , Chimonas and Hines , Testud ). 

(34) 

9. Turbulent motion in the 55 - 60 km region of the atmosphere (Murphy et al . ) 

10. In situ generation of gravity waves from non-linear interactions of tidal 

(35) 

modes (Justus '). 
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TABLE 3. POSSIBLE MECHANISMS FOR THE DISSIPATION 
OF ATMOSPHERIC GRAVITY WAVES 

1. Mode coupling i.e. , nonlinear processes in which the energy goes from one 

(35) 

wave into others (Spizzichino ) . 

2. Frictional damping* produced by viscosity and thermal conductivity Pitteway 
and Hines ^). 

(37) (27) 

3. Generation of turbulence (Hodges , Tchen ). 

(39) 

4. Vibrational relaxation (Bauer and Podney - - - see also item 7 in Table 2)* 

(39) 

5. Plasma damping effects (Liu and Yeh ). 

(39) 

6. Radiation damping (Gille ). 

7. Chemical relaxation, or the shift of the 0^ = 2 0 balance under the 

influence of the sound wave. 

8. Diffusion (in the case of a gas mexture) . 

“5/3 

9. Cascade of energy up the spectrum to larger scales because of a k two 

dimensional horizontal spectrum (see Section 5 of this report, see also 
item 1 in this table) . 
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been measured simultaneously to establish whether the observed wave-like 

structure has either the precise horizontal to vertical scale ratio for a given 

frequency wave length combination as required by the gravity wave dispersion 

relations, or the appropriate relative phase between the velocity, pressure, 

density, and temperature variations. Since the total variations are made up of 

such phenomena as planetary wave, synoptic variations, tides, gravity wave, and 

turbulence, it is not surprising that the gravity wave components lack accurate 

observational resolution, in view of the difficulties in obtaining the spatial 

(42) 

or time resolution adequate for their detection and delineation 

In view of the limitations of observation data available, the best approach 

for establishing values of parameters is probably a statistical approach such 

(43 44) 

as that applied by Woodrum, Justus, and Roper * to chemical release and 
meteor wind data. This approach considers the statistical properties of the 
non-diurnally-repeating components of wind (or temperature, etc.) and allows 
estimation of amplitudes and probability distribution. Structure functions of 
the velocity differences also indicate possible frequency spectra, and vertical 
wave length spectra. Similar structure functions of velocity differences over 
spatial separations indicate possible horizontal wave length spectra. This 
type of statistical approach is desirable since isolated single gravity wave 
components will not always be encountered. Collections of different gravity 
waves would not display the separate characteristics of the individual compon- 
ents, and the possibility of strong interaction between atmospheric tides and 
gravity waves has been suggested from both meteor and chemical release data^^*^"^ 

/ o C. \ 

and apparently confirmed by the radar meteor studies of Spizzichino v 

The purpose of this report is two fold: first to analyze the existing 

data in the 50 - 200 km range for irregular variations which could be due to 
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gravity waves, and secondly to determine methods of evaluating the effects of 
these observed irregular winds, densities, pressures, and temperatures on the 
stability, control, and dynamic heating of space craft such as the space shuttle 
vehicle. The following sections present results of measurements which establish 
the magnitude, statistics, and spatial and temporal structure of the irregular 
atmospheric variations. Section 7 of this report develops analytic methods, 
based on an analogy of the irregular motion field with axisymmetric turbulence, 
which allow measured or model correlations or structure functions to be used 
to evaluate the effective frequency spectra of scalar and vector quantities on 
a space craft moving at any speed and at any trajectory elevation angle. 
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2. DATA ANALYSIS INTERPRETATION 


Structure Functions . An alternate form of the correlation function known 
as the structure function was first used extensively by Russian meteorologists 
in the analysis of turbu^/Lnce. The structure function of a statistically stat- 
ionary time varying process f(t) is given by 

D(x) = <[f(t + t) - f(t)] 2 > (1) 

where the structure function D depends only on the time displacement x because 

of the statistical stationarity . The angle brackets in (1) denote averaging. 

Structure function analysis can be applied to wave phenomena also. Consider 

the process f(t) to be a wave of frequency to and amplitude A, i.e. f(t) = 

A sin (cot + a), where a is a phase angle. If the averaging process is taken to 

be integration over any whole number of periods, where the period T = 27r/to s 

2 2 

then the mean square value of f(t) is <f > = A /2. It is also easy to deter- 

2 

mine that the structure function in this case is given by D(x) = A (1 - cos t:x 

Thus the structure function of a single component wave field has the properties c 

2 

being proportional to x for very small time displacements x (since the leading 

2 2 

term in 1 - cos cox is to x /2) and being periodic with the same frequency co as 
the wave. Notice that the first maximum in D(x) would occur at x = tt/co = T/2 

so the first maximum occurs at the half period value. The value of the maximum 
is 2A^. 

The structure function analysis of waves can easily be extended to a set 
of waves which are part of a harmonic series. Consider f(t) to be periodic with 
period T and to be made up of N different waves each wave n having a frequency 
co^ = 2nTr/T, i.e. 
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f(t) = E A sin (u> t + a ) 
u n n 


( 2 ) 


where the summation is from n = 1 to n = N. Using the same type of averag- 

ing by integrating as on the single wave, one can determine that the mean square 
value of f is given by 


2 2 
<f > = £ A /2 

n 


(3) 


and that the structure function is given by 


D(t) = £ A (1 - cos a) x) 

n n 


(4) 


Thus again the behavior at small time displacement is proportional to t , since 

2 2 2 

the leading term in (4) is £ x /2, The structure function in this case 

is again periodic with period T and has components with amplitudes related to 
the square of the amplitudes of the corresponding frequencies in the process 
f(t). 

If the process f(t) is made up of sufficient components that it must be 
considered as having a continuous spectrum of amplitudes A(w) (as is the case 
with turbulence and may also be the case when gravity wave modes have signifi- 
cant nonlinear modal interaction) then f(t) can be represented by 


f(t) 


-J 


A (o)) e lt0t da) 


(5) 


(For a more rigorous discussion of spectral representation see Lumley and Pan- 
ofsky^^ and Lin^^). The mean square value of f(t) is given by 


<f > = 2 

— cqJ 


4>(oi) do) (6) 

where c}> (tn) is the spectral density of the mean square variations of f(t). The 
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structure function is given by 


D(x) 



(1 - cos u>t) da) 


( 7 ) 


Again the initial variation of D(x) for small time displacements is proportional 

_ 2 . 
to t since 


D(t) ■ [ 4> (oj) dm] x 0 (8) 

o^ 

However, at large time displacements D(t) is not periodic but approaches a con- 

2 

stant value of 2<f >. This is seen readily from the fact that D(t) is related 
to the correlation function p(x) by 


D(t) = 2<f 2 > [1 - p(x)] (9) 

for all stationary processes f(t). Hence, because p(x) approaches 0 for large 
x, then 

D(x) = 2<f 2 > x (10) 


In addition, if there is an appreciable range of frequencies over which <j) (w) 
varies as a power law 

<Hw) = c o)' a (li) 

then there is a range of time displacements over which 

r 00 

D(x) Z C[ x a (1 - cos x) dx] x a 1 (12, 

O'* 

where x = ojx and the integral in brackets is some definite constant value. 

Thus observed power law behavior of the structure function (i.e. a - 1) can be 
related to the power law of the spectral density (i.e. -a). 
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Daily Difference Analysis . A technique for resolving irregular variations 

in atmospheric winds or other parameters has been previously developed by 

(43 44) 

Woodrum and Justus * for application to chemical release data. This tech- 
nique does not rely on explicit extraction of tides and prevailing values and 
does not require large amounts of data. Hence it is well suited for the study 
of irregular variations at high altitudes where data is limited in quantity. 

As an example of the application of this technique consider two vertical 
profiles F^(z) and F 2 (z) of a parameter F(z) where F might be a wind component, 
pressure, density, or temperature. Suppose that the two vertical profiles were 
measured at the same time of day and within a few days of each other. The 
assumption is made that the two profiles are represented at all measured heights 
z by the relations 


^(z) = F Q (z) + fj_(z) 

F 2 (z) = F o (z) 4- f 2 (z) (13) 

where F (z) is a profile made up of the prevailing value, tidal variation, and 
any other diurnally repeating component of F. The profiles f^(z) and f 2 ( z ) are 
the profiles made up of the gravity wave components and any other non-diurnally 
repeating components such as planetary waves, synoptic variations, and varia- 
tions of the tidal parameters from day-to-day. If the original data measuring 
technique did not filter out the turbulent components these would also appear 
as a contribution to f^(z) and £ 2 ( 2 ), and if the number of days interval between 
the measured profiles becomes large then seasonal trends will also make a contri- 
bution. It is readily apparent from (13) that the difference in the measured 
values at any height z is equal to the difference in the irregular values at 
that height, i.e. 
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AF(z) 


F 2 (z) - F 1 (z) 


f 2 (z) - £^(z) 


(14) 


Two further assumptions are now made: first that the irregular variations are 

statistically stationary and vertically homogeneous, and secondly that the irre- 
gular variations at the two measurement times are uncorrelated with each other. 
The second assumption is valid only if the irregular variations are made up of 
components whose periods are small compared to the interval between the measure- 
ments (e.g. gravity waves of up to a few hours period). If longer period phe- 
nomena, say planetary waves, of a predominate period of n days make a signifi- 
cant contribution then they should cause peaks in the average AF's, at multiples 

of n/2 days. In the present analysis the AF's are considered stationary. The 

(44) 

analysis of chemical release winds did not show significant variations of 

average AF's with the number of days difference, and the study of this question 
for variations at lower altitudes will be the subject of later investigation. 

Under the assumptions outlined above it can be shown that certain proper- 
ties of the data difference values can be related to the irregular variations. 

For example consider an average (at the same altitude) over an "ensemble 11 of 
several profile pairs. The mean square difference is thus given by 

< [ AF(z) ] 2 > = <f 2 2(z)> + <f i 2 ( z ) > - 2<f x ( z ) * 2 (z)> = 2<f 2 (z) > (15) 

where the last step is possible because of the independence and stationarity of 
f^(z) and f 2 (z). The development leading to (15) is also taken as justification 
that the probability distribution of the data differences (divided by /~2) is 
the same as the probability distribution of the irregular variations. However, 
this is strictly true only if the distribution is Gaussian. 

The vertical structure function of the data differences is also related 
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to the structure function of the irregular variations. The vertical structure 
function of the data differences is given by 

D af (?) = <[AF(z + C) - AF(z)] 2 > (16) 

where, again, the angle brackets represent an "ensemble" average. Substitution 
of (14) and rearrangement produces 

D af (t) = <[f 2 (z + c) - f 2 (z)] 2 > + <[f x (z + 0 - f x (0] 2 > 

-2<[f 2 (z + C) - f 2 (z)][f 1 (z + O - f 1 (z)]> (17) 

Application of the assumptions of independence and vertical homogeneity of the 
f T s results in 

D af (t) = 2<[f(z + C) " f(z)] 2 > = 2D f (c) (18) 

In other words, the vertical structure function of the data differences is 
twice the vertical structure function of the irregular variations. 

In addition to vertical structure functions , horizontal structure functions 
were also calculated by differencing data from the up and down trail trajector- 
ies of rockets or from closely spaced simultaneous observations. It should be 
noted that the assumption of independence of the irregular variation profiles 
would not be valid in this case, however, and this type of difference therefore 
yields a structure function which has significant variation with the distance 
between the places of evaluation of the two profiles. However, since the 
effects of tides and other phenomena are not removed, the horizontal structure 
functions yield valid results on truly irregular variations only if the dis- 
tance between profiles is kept small compared to the horizontal wave lengths 
of these tides and other phenomena. It should also be noted that the hori- 


ft 
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zontal structure function used for velocity is actually the structure function 
of the horizontal velocity magnitude 

D(r) = <[u(x + r) - u(x)] 2 > + <[v(x + r) - v(x)] 2 > (19) 

where u is the eastward velocity component, v is the northward component, x is 
the horizontal vector location of one velocity evaluation, r is the vector dis- 
placement between the locations of velocity evaluation, and r is the magnitude 
of the vector r. The definition (19) is dependent on the magnitude r and not 
the direction of the displacement and no separate considerations of longitudinal 

and transverse structure functions are required. (A longitudinal structure 

2 2 

function would be of the form <[u (x + Ax) - u(x)] > and a transverse structure 

2 2 

function would be of the form <[u ( y + Ay) - u(y)] >. 

Harmonic Analysis . In addition to a large irregular component, the winds 
and thermodynamic variables in the upper atmosphere are composed of regular 
variations produced by diurnal, semidiurnal, and other tidal components, and 
possibly long-period variations such as planetary waves. At the present time 
there is no method of resolving the values determined from a single profile 
into these various components. However, a series of profiles determined at a 
sequence of closely spaced times can be used to resolve the various components. 

A sequence of measured values at a given altitude are fit by a least squares 
process with the function 

F(t) = Aq + A ^2 sin(27rt/12 + <j> 12 ) + A, ^ sin(27rt/24 + <j> 2 ^) (20) 

where the time t is measured in hours, Aq is the mean value of the parameter, 

A 12 and are the amplitudes of the 12 hour, and 24 hour period tides, and 
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and <}>24 are phase angles for those tidal components. The A’s and cf> T s are 
the parameters determined by the least squares process (indirectly by first 
evaluating coefficients of sine and cosine terms equivalent to (20)). 

The computed tidal parameters were evaluated at each kilometer altitude 
point and were found to form smooth profiles of variation with altitude. After 
the mean values and tidal components have been determined from a series of pro- 
files, then these values may be subtracted from the observed values to yield 
residual components. These residual values can be evaluated at each altitude 
for each of the profiles in the sequence and profiles of irregular variation 
values constructed. Each of these irregular parameter profiles can be used to 
evaluate vertical structure functions and values from different profiles in the 
time sequence can be differenced to computed time structure functions. 
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3. PROBABILITY DISTRIBUTIONS AND MAGNITUDES 


Distribution Statistics , The probability distribution p(f) of a parameter 

f(x, t) is defined so that p(f Q ) df is the probability of observing a value of 

the parameter with a value between f and f + df. The usual normalization is 

o o 

to unit probability for observing a value of f in the range to + This is 
expressed mathematically by the relation 


p(f) df - 1 


The average value of f is found by the equation 


f = <f> 


f p(f) df 


and the probability distribution function can be characterized by the moments 
y, about the average. The moments are given by 


P k = <(f - f) k > = [ (f - f) k p(f) df 


Of particular interest are the standard deviation a, the skewness a, and the 
kurtosis 3 which are defined by 


O = \lr 


a = u 3 /a~ 


' p 4 /a 4 


The Gaussian probability distribution is given by 


p(f) - exp[-(f - f) 2 /2a 2 ] / (271 ) 1/2 * 


and is thus completely characterized by its mean value and standard deviation. 
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The skewness (as with all odd moments) for the Gaussian distribution is zero. 

The kurtosis for the Gaussian distribution is 3. 

Characteristics of the probability distribution of irregular winds at chem- 

(44) 

ical release altitudes have been reported by Justus . For this study data 

on winds, pressures, densities, and temperatures were gathered from a large 

number of sources and analyzed in the same fashion, (i.e. as outlined in Section 

2). These included the Meteorological Rocket Network (MRN) data for the years 

1964 - 1969^^; NASA Goddard grenade and pitot tube data for 1960 - 1969^^; 

chemical release winds from Eglin AFB, Florida 1962 - 1971^^, from Barbados 

and Yuma^^, Wallops and other locations . ROBIN falling sphere data^"^ * and 

( 54 } 

other rocket and falling sphere data v ' . In addition to these collected data, 
individual data sets from many sources were added to the computed probability 
distributions^^»20,54 ^1) . These data include satellite measurements, meteor 
winds, and ion probe observations. 

Figures 1 through 5 show some representative results of the probability dis- 
tribution analysis. The ordinates of these graphs are the number of observed 
values. Bars forming the shaded area represent the observed values and the large 
dots are the expected number based on zero mean and the observed standard devia- 
tion in a Gaussian distribution. Relative probabilities may be obtained by 
dividing the number of observations by the total number observed, given on each 
graph. The abcissa on the velocity plots. Figure 1 and 2, are the equivalent 
irregular wind magnitudes determined by the observed data value velocity differ- 
ences divided by the square root of two (cf. equation (15)). For pressure, 
density, and temperature distributions the abcissas are relative fluctuation 
magnitudes i.e. p , /p Q , p f /p Q , T f /T q where the fluctuation values are determined 
by the observed data differences divided by v^2 


.1 
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(26) 


P* = (P 2 " P 1 )/» / 2 _ 

p’ = (p 2 - p 

T' - (T 2 - T X )//F 

and the mean values are determined by the average of the observed data values 

p o = (p l + p 2 )/2 

P 0 ° (pj + P 2 )/ 2 (27) 

T q = (T x + T 2 )/2 

The observed magnitudes, as characterized by the standard deviations of 
the distributions, showed latitude and altitude variations which are discussed 
later in this section. The computed values of skewness parameter a showed no 
variation with height or latitude. Values computed for the kurtosis g did 
show a significant altitude variation discussed below. Table 4 shows a summary 
of the calculated a values for the irregular winds, pressure, density, and tem- 
perature. The values of a, a, and g were computed from the observed distri- 
bution in histogram form, with Sheppard T s corrections applied, rather than 
from each individual observed value. In order to test the significance of 

the deviations of a from the expected Gaussian value of zero the test explained 
f 62 ^ 

by Keeping^ ' was employed. The value used in the test for the number of data 
was the estimated number of independent data values rather than the total num- 
ber observed. The number of independent values is less than the total for two 
reasons. First the data values at adjacent 1 km height intervals are not in- 
dependent but are correlated with each other, and secondly if a set of data 
from n consecutive days are to be differenced then n(n - l)/2 differences can 
be formed, but only n - 1 of these differences are independent of each other. 


i 
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Table 4. Average Skewness a of 
the Probability Distribution of the 
Irregular Variation 


Irregular 

Variation 

Parameter 

Total 
Number of 
Observations 

Number of 

Independent 

Observations 

a 

p' 

4175 

1252 

-0.11 

p' 

4151 

1252 

-0.09 

T' 

4413 

1432 

-0.22* 

u 

(eastward) 

9559 

2359 

+0.13* 

V 

(northward) 

9078 

2283 

-0.10* 


Values which are significantly different (at the 1% level) 
from the zero value expected for a Gaussian distribution. 
Significance tests are those given by Keeping^ 'and based 
on the number of independent observations. 
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To account for the correlation over height it was assumed that data values 
separated by 5 km of altitude could be taken as independent for purposes of 
computing the number of independent observations. 

Three of the skewness values in Table 4 were computed to be significantly 
different from zero. A positive a indicates an excess of large positive 
values and a deficit of large negative values compared to a symmetric distri- 
bution. The reverse of this would be true for a negative a. The largest 
magnitude skewness is for the irregular temperature distribution. Because of 
the small a of the temperature there were few intervals into which the values 
fell and the histogram technique, even with Sheppard f s corrections, probably 
does not accurately estimate the higher moments. Of the two remaining a f s 
which are significantly different from zero, one is positive and the other 
negative. Therefore, it is not clear what if any significance should in fact 
be placed on these skewness results. 

The values computed for the kurtosis g showed a significant altitude 
variation, as reflected in Figure 6. The dotted area represents the approxi- 
mate range within which 3 would not be significantly different from the expected 
Gaussian value of 3. The notation used in Figure 6 is u for the eastward wind 
component and v for the northward wind component. The values are consistently 
high-to-normal below 100 km and low-to-normal above 100 km. The extremely 
large values found for the irregular temperature may be due to the same prob- 
lem of low o and small number of intervals, which also led to high magnitudes 
for the skewness. 

High values of 3 indicate a distribution with a high central peak and 
broad "skirts”. A low 3 indicates the reverse of this. The height variation 

(44) 

of 3 is consistent with the previous measurements by Justus of low g, for 
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the irregular winds, because they were for the region near 100 km where low 3 
values are indicated. Also shown on Figure 6 are the 3 values of the distri- 
butions of irregular winds determined as residual winds from meteor data from 
Garchy, France and Durham, New Hampshire^^ after resolution of the tidal 
winds by a data reduction method developed by Groves v . The meteor data 
indicate near Gaussian values for 3 in the 80 - 100 km region, in slight dis- 
agreement with the other data, but in good agreement with Russian meteor results 
computed by Pokrovskiy, et al.^^. 

The high 3 values of the MRN data at 50 - 60 km persisted even when the 
data was devided by site location. The 3 values computed for each site separ- 
ately and then averaged were approximately the same as those plotted in Figure 
6, except for the relative temperature fluctuations. The site averaged 3 of 
the T f /T distribution was 6.06 at 50 km and 9.34 at 60 km, somewhat lower than 
the values of Figure 6. 

Magnitude Variation with Height . The magnitude of the irregular variations, 
as represented by the standard deviation a of the probability distribution, 
showed significant height variation. The measured height variation of the 
irregular winds is seen in Figure 7. The data which went into Figure 7 included 
the 1969 MRN data from 8 sites: Fort Sherman, Cape Kennedy, Ascension Island, 

White Sands, Fort Greeley, Point Mugu, Barking Sands, and Fort Churchill, plus 
all of the data from other sources. The magnitude determined by the daily 
difference method from the MRN and other data are shown by the data points 
connected by solid and dashed curves. In Figure 7, as elsewhere, u is the 
eastward wind component and v is the northward component. At 50 and 60 km the 
irregular winds are definitely anisotropic, with the eastward component domi- 
nating. At 70 km and above no significant anisotropy exists. In view of the 
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Figure 7. Height variation of the magnitude of the irregular winds. u is the 
eastward wind component and v is the northward component. 





size of the error bars above 120 km, the apparent dominance of the eastward 
component at those heights is not considered to be significant. The increase 
of irregular wind magnitude with height is rather sharp between 90 and 100 km 
as measured by the daily difference data. A rapid though not quire so sharp 
increase at these altitudes is also seen from the Durham meteor data irregular 
winds. Although the Garchy data agrees in general magnitude, the sharp rate of 
increase near 90 km does not appear in the Garchy meteor data. A still differ- 
ent value near 90 km was obtained from Russian meteor winds from Kazan ^4) _ 

The Kazan value appears to be more representative of the higher altitude data 
than the Durham or Garchy data. These variations lead one to conclude that the 
height variation curve of Figure 7 given by the MRN and other data may be valid 
in detail only for the North American continent. 

Linearized gravity wave theory with no dissipation predicts that the 

2 2 

gravity wave kinetic energy density p Q <u + v > remains constant with height. 

It is obvious from Figure 7 that considerable dissipation takes place above 
100 km since the irregular wind magnitude remains relatively constant while 
the atmospheric density decreases rapidly. However, significant dissipation 

also occurs between 50 and 100 km despite the increasing irregular wind magni- 

—4 

tude. For example, the density decreases by a factor of about 5 x 10 from 
50 to 100 km while the mean square irregular velocity increases by a factor 
of only about 8. 

The height variation of the irregular pressure, density, and temperature 
is shown in Figure 8. The data are expressed as variations relative to the 
mean atmospheric parameters. The declining values above 110 km confirm the 
fact that significant dissipation of the irregular variations occurs above 
about 110 km. The dissipation above 100 km is not unexpected in view of 
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RMS Value, % 

I 

Figure 8. Height variation of the magnitude of the relative pressure p / p Q , 
density p'/P 0 > and temperature T ' /T^ . 
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earlier theoretical work on the dissipation of gravity waves^*^^ and experi- 

(43 44 ) 

mental observations of the height variation of irregular winds v * ' in this 
region. However the more-or-less continuous dissipation from 50 to 100 km is 
somewhat unexpected. 

Figure 8 shows that the relative temperature variation magnitude between 
50 and 70 km is significantly smaller than for the relative pressure or density, 
which are about equal. Between 80 and 100 km the relative pressure and’ temper- 
ature variation magnitudes becomes about equal and smaller than the relative 
density variations. Above 100 km the relative density variations become smaller 
than those for the relative temperature. 

At first one might think that the Boussinesq approximation should apply to 
the pressure, density, and temperature variations and hence that the relative 
pressure variations should be much smaller than the relative density and tem- 
perature. However, as shown by Dutton and Fichtl^^, for motions of vertical 
scale not small compared to the scale height, the magnitudes of all three 
parameters are comparable. It will be shown later that the vertical scales are 
comparable with the scale height. Therefore, the relative values of the irre- 
gular thermodynamic parameters are related by 

p'/p o = p'/p Q + T'/T o (28) 

and in the mean square they would be related by 

<(p’/p o ) 2 > = <(p’/p o ) 2 > + <t’/t o ) 2 > 

+ 2<(p*/p o )(T 7T o )> (29) 

where the last term is basically the correlation between the density and tem- 
perature variations, which in the case of a single component wave train would 
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Figure 9. Comparison of observed height variation of p /p and T /T with 

(67) T ® ^ 

Miller T s V ' theoretical curve of p /p^ based on measured irregular 

wind magnitudes. 




be dependent on the relative phase between the density and temperature wave 
variation. The changes discussed above in the relative magnitudes of the 
thermodynamic variables of Figure 8 therefore indicate a changing phase 
relationship with height. This changing phase relationship could be due to non- 
linear intermodal interactions or to selective model absorption of a complex 
gravity wave field. In either case it would be difficult to justify the chang- 
ing relative magnitudes with the assumption of a simple mode gravity wave field. 

A comparison of the height variation of the irregular winds and irregular 
density shown in Figure 9 gives good agreement with a theory developed by 
Miller^^, which predicts relative density variations proportional to rms ve- 
locity variation by a factor which varies with height . The theoretical values 
of P ! /p o are related to the mean square irregular winds by the formula 

IpVpJ = {[(y - D/c 2 + (l/g T 0 ) dT o /dz] <u 2 + V 2 >} 1/2 (30) 

The U.S. Standard Atmosphere of 1962 was used to evaluate the speed of sound c, 
the acceleration of gravity g, the temperature T q and its height gradient 
dT Q /dz. The ratio of specific heats y was assumed to be constant at a value of 
1.4. The agreement between the observed p f /p o variation and that predicted by 
the theoretical formula is remarkably good, considering the uncertainties intro- 
duced by the use of standard atmosphere values. However, it does appear from 
the results shown in Figure 9 that the theoretical formula corresponds more 
nearly with the magnitude of the irregular temperature T f /T Q than the density 

p’/p 0 . 

Magnitude Variation with Latitude . Only the Meteorological Rocket Network 
data were available in sufficient quantity to study latitudinal variations. 
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The data which were used are 1964 - 1969 results from the MRN sites listed in 
Table 5, The average magnitudes, as determined from the daily difference method 
standard deviation of the distribution, are plotted in Figure 10 as a function 
of the observing station latitude. The curves identified as P and T are rel- 
ative pressure p 1 /p Q and relative temperature T f /T Q in percent. The latitude 
variation of the relative density was essentially identical with that for relat- 
ive pressure. The irregular velocity magnitudes are identified as U for east- 
ward component and V for the northward component. 

At both altitude levels in Figure 10 the relative temperature variation 
magnitudes are considerably smaller than the relative pressure or density. This 
is in keeping with the results previously shown in Figure 7. The magnitudes of 
the relative temperature variations are so small (1 - 2%) that one at first 
questions the validity of the interpretation that these represent atmospheric 
variations and not just observational error. First it should be pointed out 
that the daily difference method involves observations at the same time of day 

/£Q\ 

and therefore extraneous variation due to radiation errors ^ from measurements 
at different solar elevations cannot be a source of error. Secondly results of 
the study of rocketsonde repeatability indicate temperature measurement 
accuracy of about 1°C which at 50 - 60 km is about 0.4% (velocity accuracy is 
estimated at about 3 m/s) . Thus the observed relative temperature variation 
magnitudes of Figure 10 are between 2.5 and 5 times larger than the observational 
limit. Thirdly, the trend of increasing magnitude of variation with increasing 
latitude is consistent with the observed increased variability of the atmosphere 
at high latitudes w hile it would be difficult to explain an increasing obser- 
vational error with increasing latitude. 

The general trend of increase with latitude is seen in all of the varia- 
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TABLE 5. MEN SITE LOCATIONS 


Site Name 

Ascension Island, A.F.B. 
Fort Sherman, Canal Zone 
Barking Sands, Hawaii 
Cape Kennedy , Florida 
White Sands, New Mexico 
Point Mugu, California 
Fort Churchill, Canada 
Fort Greely, Alaska 


Latitude Longitude 


7° 

59' 

s 

14°' 

25' 

W 

9° 

20 ' 

N 

79° 

59' 

W 

22 ° 

02 ’ 

N 

159° 

47' 

W 

ho 

00 

o 

27’ 

N 

o 

o 

00 

32' 

W 

32° 

23' 

N 

106° 

29' 

W 

34° 

23' 

N 

119° 

07' 

W 

58° 

44 < 

N 

93° 

49' 

W 

64° 

00 ' 

N 

145° 

44' 

W 
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bles plotted in Figure 10. Two Exceptions to the trend are the sharp peak at 
32.5° N latitude (the White Sands site), and the sharp drop in the 60 km data 
at 64° N latitude (the Fort Greely site). In an attempt to discover a possible 
cause for the peak value from the White Sands results a possible geographical 
cause was noted in the topography west of White Sands. Figure 11 shows several 
southwestern states and the upper portion of Mexico. The location of White 
Sands is labeled A. The shaded area shows the approximate region which is 
fairly densely populated with mountains with peak height of approximately 2 km 
or higher. The western portion of this area is made up of the Sierra Nevada 
Mountain Range, the middle portion is the Colorado Plateau Area and the Eastern 
section is a portion of the Rocky Mountains. Some peaks in the Sierra Nevada 
Range are identified in Figure 11: B is Mt. Whitney (4418 m) , C is Toro Peak 

(2653 m) , D is Mt. Palomar (1867 m) , E is Santa Catalina (3089 m) . Immediately 
East of the Sierra Nevadas is an area which is at or below sea level: F iden- 
tifies Death Valley (-86 m) , G is the Salton Sea (-72 m) , and H is the Gulf of 
California. Any wind blowing toward White Sands from the WNW, W, or WSW would 
encounter this mountain-followed-by-lowlands terrain. Since White Sands is 
approximately 1000 km from the Sierra Nevada Range, then the dominant period 
of waves launched by those mountains reaching 50 - 60 km height at the location 
of White Sands should have periods near 150 min. or longer i.e. slightly more 
than two hours. This result was calculated using the method employed by 
Hines on gravity waves generated by nuclear explosions. This value is in 

general agreement with the periods of waves inferred later in ths report from 
time structure function analysis. The reason why the Fort Greely values, esp- 
ecially at 60 km, are low is not clear at this time. However, the results of 
(72) 

Theon et al . ' have shown that wave structure evident in the 40 - 90 km 
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region over Pt. Barrow (71° N) in the winter (January - February) is entirely 
absent in the summer (June - August) . This result indicates strong seasonal 
variation in the magnitude of upper atmospheric waves at high latitudes, possi- 
bly due to variation in the propagation characteristics of the atmosphere between 
the ground level and 40 km. 
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4. VERTICAL STRUCTURE FUNCTION DATA 


The vertical structure function is defined as 

D(C) = < [V (z) - V (z + ?] 2 > (31) 

where V may be velocity component, pressure, temperature or density, z is the 
height coordinate and ? is the height displacement. 

Data . Profiles of residual winds, pressure, temperature and density were 
obtained from Meteorological Rocket Network Data by differencing pairs of wind, 
pressure, temperature, and density profiles according to the method discussed 
in section 2 (see equation (15)). The pair members were chosen such that both 
profiles occurred within ten minutes of the same local time of day and not more 
than fifteen days apart. This procedure was performed on data from the eight 
locations listed in Table 5. The time coverage of the data ranged from 1964 to 
1971 and the height range was from 45 km to about 65 or 70 km. 

In addition, similar profiles or residual winds were obtained from chemi- 
cal release wind data taken at Eglin AFB in Florida , GCA data^“^ and gun 

probe data from Barbados and Yuma^^ . The time coverage of this data was from 
1959 to 1971 and the height range was from about 80 km to 140 km. The chemical 

release data also includes some vertical winds which were used in the vertical 

i - (50,73,74) 

structure analysis 

The ROBIN falling-sphere data taken at Eglin AFB on May 9-10, 1961, was 
(53) 

used also . The height range of this data was from about 45 to 65 km. 

Irregular winds were also obtained from the ROBIN data by subtracting the 
tidal components determined by harmonic analysis from the raw data. 

Analysis. The vertical structure functions of the residual winds were 
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computed for each residual wind profile by taking differences between values 
from the profile at altitudes different by. the displacement amount and averag- 
ing over the profile. Similar analysis was performed for the pressure, temper- 
ature, and density data. 

The MRN data were divided into groups according to the site location from 
which the data were taken, and data from each site was analyzed separately. 
Representative profiles of the vertical structure functions of these data are 
shown in Figures 12 and 13. At first the residual profiles from a particular 
site were divided into two height groups of 45 to 55 km and 55 to 65 km, and a 
structure function was calculated for each height group. However, the struct- 
ure functions for the two groups were not significantly different; and thus, 
the height groups were recombined into one group. 

Figure 12 shows the structure function for two profiles of residual north- 
south winds and two profiles of residual east-west winds. All four functions 
show the same general characteristics. For small vertical displacements the 
structure function increases approximately linearly on log-log scale and then 
levels off. A power law of the form D(£) = const. (£) n was fitted to each of 

the vertical structure functions of the MRN data over small vertical displace- 
ments £. The averages of the corresponding exponents n were 1.06 for the 
residual north-south winds and 0.94 for the east-west winds. 

The average leveling-off values of the structure functions were found to 
2 2 

be 197 m /sec for the residual north-south winds, corresponding to a Z 10 
2 2 

m/sec, and 237 m /sec for the residual east-west winds, corresponding to 
Z 11 m/sec, in reasonable agreement with the values at 50 and 60 km in 
Figure 7. 

Also, characteristic with most of the vertical structure functions was 


41 


<Av >, (•"/• 



Figure 12. Vertical structure functions of MRN irregular wind data determined 
by the daily difference method. 
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the additional peak for a vertical displacement of about 18 km. This peak could 
be caused by a wave motion with a dominant vertical wavelength of about 36 km. 

The amplitude of this motion seems to . vary from about 5 to 10 m/sec. 

The horizontal arrow in the graph of the structure function for the north- 
ward winds at the Cape Kennedy site gives an example of the position of the level- 
2 

ing-off point (2a ) which was calculated from the daily-difference analysis. 

Actually, the daily-difference analysis calculates an average for the leveling- 

2 

off value and would include the peak in the average. Thus, the value of 2a 
calculated from the daily-difference program should be slightly larger than the 
value calculated from the vertical structure functions as is shown. 

Figure 13 shows the vertical structure functions for two profiles of rel- 
ative temperature deviation and one profile each of relative pressure and den- 
sity deviation. The function for the temperature seems to behave similarly 
to the structure functions for the wind components in Figure 12. However, the 
vertical structure functions for the pressure and density have entirely differ- 
ent detailed characteristics. They do seem to have a linear log-log increase 
for small vertical displacements, but they have no leveling-off point. Some 
of the functions seem to peak at vertical displacements of 18 or 19 km similar 
to the peaks observed in the velocity structure functions whereas others con- 
tinue to increase out to 20 km, the largest displacement calculated. 

A power law D(c) = const. (c) n was fitted to the vertical structure 
functions for the temperature, pressure, and density data over small vertical 
displacements. The averages of the exponent n are 0.82, 1.09,. and 0.65 for 

the pressure, temperature, and density respectively. The average leveling- 

2 

off value of the temperature is 5.0(%) corresponding to a = 1.6%, again in 

good agreement with previous results in Figure 8. 
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Figure 14 shows the vertical structure functions for the chemical release 

wind data. These data were broken into two height intervals, 80 to 110 km and 

110 to 140 km. The functions show a linear log-log increase with increasing 

vertical displacement and then leveling-off. The leveling-off values are much 

larger than those for the MRN data as one would expect due to the increased 

2 2 

altitude. In this data the leveling-off value is about 3500 m /sec corres- 
ponding to a Z 42 m/sec, in reasonable agreement with Figure 7 for this height 
range. 

Figure 15 shows the vertical structure functions for the vertical winds. 

The raw vertical wind data were used in the computations with no daily differ- 
encing. However, since essentially all of the vertical wind should be the irre- 
gular component this is considered appropriate. Both functions, for data below 
and above 110 km, show a linear log-log increase over small vertical displace- 
ments with slopes of 1.32 and 1.52 respectively. The function for the data 

above 110 km does not have a definite leveling-off point, whereas, the function 

2 2 

for the data below 110 km levels off at a value of 250 m /sec corresponding to 
a 2 11 m/sec. This a is much smaller than the above calculated a- value for 

the horizontal winds in the same altitude range. Again, this is as expected 
because the vertical amplitudes of atmospheric gravity wave motions are expected 
to be smaller than the corresponding horizontal amplitudes. 

The ROBIN daily difference data and ROBIN irregular wind data from harmonic 
analysis were analyzed similarly and the results are shown in Table 6. The re- 
sults of the ROBIN data were not averaged with the results of the MRN data be- 
cause the resultant values were significantly different. It is suspected that 
something peculiar is occurring in the ROBIN data and consequently, the ROBIN 
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values were kept separate. 

Table 6 also gives a summary of the results of the vertical structure 

function analysis and other results. 

(44) 

Justus analyzed chemical release data from Yuma to obtain a vertical 
structure function via an analysis method which is different from the method used 
in this report. His results are shown in Figure 16 and they agree with the a 
results summarized in Table 6, but the exponent n = 1.6 is slightly larger. 

If the exponent values n of Table 6 are interpreted via equation (12) as result- 
ing from a continuous power law spectrum, then spectra of the vertical wave num- 
ber k would follow a power law <f> (k) « k ^ + n \ The tendency of the Table 6 
values of n to increase with altitude would thus indicate a changing spectrum 
with height. 

The dominant vertical scale of the irregular motion can also be observed 
from the vertical structure functions. The half wavelength of the dominant 
vertical scale is given as the minimum height displacement for which the verti- 
cal structure function begins to level off. An average value for the dominant 
vertical scale of the MRN data (data heights from 45 to 65 km) is about 7 or 
8 km. The same average value was found also for the chemical release data which 

is at a height of around 110 km. These results are in agreement with theoretical 

(27) 

predictions by Tchen , who stated that the minimum vertical scale of internal 
gravity waves should be proportional to the pressure scale height. The scale 
height for the height regions of concern is fairly constant at about 8 km. 

Earlier experimental support of the correlation between scale height and mini- 
mum vertical scale was given by Justus in his analysis of chemical release 
data. However, the structure function scale of 7 km at the 45 to 65 km level 
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TABLE 6. SUMMARY OF STRUCTURE FUNCTION INITIAL 
EXPONENTS n AND MAGNITUDES a INFERRED FROM LEVELING-OFF VALUES 


Data 

Height 

Range 

(km) 

V 

ns 


V 

ew 


Temperature Pressure Density 

V 

z 



n 

a 

(m/sec) 

n 

a 

(m/sec) 

n 

a 

n 

n n 

(m/sec) 

MEN 

45-65 

1.06 

10 

0.94 

11 

1.09 

1.6 

0.82 

0.65 


ROBIN (Daily Diff.) 

45-65 

0.69 

18 

0.60 

18 

0.56 

— 

1.32 

0.72 


ROBIN (Harmonic 
Analysis) 

45-65 

0.84 

— 

0.68 

— 

1.52 

6 

1.48 

1.40 


Chemical Release 

80-110 

1.32 

42 

1.32 

47 




1.32 

11 

Chemical Release 

110-140 

1.08 



1.44 

42 




1.52 





10 
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Figure 16. Vertical structure function of chemical release data previously 


analyzed by Justus 


( 44 ) 


using a deviation-from-the-mean technique 




(-7 £ \ 

is contrary to the estimates by Webb v J and others of a vertical scale of 
about 2 km in that height range. These lower estimates have come from sub- 
jective estimates of layer thickness in the wind profiles or correlation scales 
of deviations from 5 km vertical means. The first of these may be deficient 
because of tendency for the eye to ignore larger scale undulations for the 
smaller scale "wiggles" in the profile, and the second method a priori excludes 
the computation of a scale larger than 5 km. 
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5. HORIZONTAL STRUCTURE FUNCTION DATA 


The horizontal structure function is defined as 

D(r) = <[V(x) - V(x + r) ] 2 > (32) 

where V may be the horizontal velocity (made up of the two horizontal velocity 
components), the pressure, temperature, or density; x is the horizontal coordi- 
nate vector; and r is the horizontal displacement vector. The structure 
function depends only on r, the magnitude of the vector r. 

Data . The horizontal structure function of the irregular winds can be 
estimated by taking differences in winds determined from up and down leg trails 
or by taking differences in winds determined from two simultaneous measurements 
separated by a small horizontal distance. The horizontal separation of the 
measurements provides the horizontal displacement in the structure function 
computation. Data meeting these qualifications were found in the chemical 
release data obtained from rockets fired at Eglin AFB, Florida and in the 
chemical release data obtained from the gun firings at Barbados, West Indies, 

and Yuma, Arizona^'*' * . The altitude range of these data is about 80 to 140 km. 

(53) 

In the lower altitudes, an attempt was made to use the ROBIN data , but 
there were not enough data available to obtain statistically meaningful infor- 
mation. Also, the same was true for pressure, temperature and density data. 

(78 79) 

Although some data were available^ * , there were not enough data to deter- 

mine any trends of the horizontal structure function for these quantities. 

However, the results computed here and averaged with structure function values 

reconstructed from correlation values reported by Mahoney and Boer^^, Cole 
(81) (82) 

and Kantor'' 1 , and Lettau'' did produce reasonable results. 


52 



Analysis . Figure 17 shows the horizontal structure function of the hori- 
zontal winds. Horizontal structure function data for the vertical winds was 
also obtained from the chemical release data and is shown in Figure 18. Each 
data point on these figures represents the average of several data' points which 
were grouped according to small horizontal separation intervals. The data 
point at the largest separation (450 km) represents the results of differenc- 
ing one pair of chemical release wind profiles reported by Blamont v ' . The 
line through the data points of Figure 17 represents variation according to the 
power law 

D(r) = const, (r)^^ (33) 

and is seen to be in reasonable agreement with the observed data. This is also 

in agreement with previous results reported by Justus and co-workers . 

Such a power law would correspond to a spectrum of horizontal wave numbers k 

-5/3 

of the form <p (k) * k . One of the theoretical results of two dimensional 
turbulence^^ is that a -5/3 spectrum would lead to an energy cascade up the 
spectrum to higher scales. Therefore it may be that the gravity wave field 
in the upper atmosphere, in addition to losing energy by damping and the gen- 
eration of turbulence, may also serve as a partial energy source for larger 
scale motions, e.g. planetary waves and synoptic scale variations. The re- 
sults of Figure 17 are compatible with a magnitude (a) of about 40 m/s and a 
horizontal scale (half wave length) on the order of 500 km. 

Conclusions from the vertical winds structure function in Figure 18 are 
less certain. One interpretation as shown would yield a power law of 1.06 at 
small scales and indicate a horizontal scale of about 20 - 30 km. The fact 
that the horizontal scale of the vertical winds is so much smaller than the 
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Figure 18. Horizontal structure function of the vertical chemical release wind 




scale for the horizontal winds indicates that the small scale waves have larger 
contributions from vertical components than do the large scale waves. 

The horizontal structure function of the horizontal winds in the 45 - 65 km 
range are shown in Figure 19. The three data points at the largest scales are 
from Cole and Kantor^ 81 ^ and other data points were averaged from the other 
available data. The 2/3 power law fits the data reasonably well at small scales 
and the structure function at larger scales indicates a horizontal scale of about 
100 km and a magnitude of about 10 m/s, somewhat lower than inferred earlier in 
this report. 
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6. TIME STRUCTURE FUNCTION DATA 

The time structure function is defined as 

D(At) = <[V(t) - V(t + At)] 2 > (34) 

where V may be a velocity component, pressure, temperature, or density, t is 
the time and At is the time separation. 

Data . Profiles or irregular winds, pressure, temperature and density were 
obtained by subtracting tidal contributions from raw data. The tidal contri- 
butions were found by performing a harmonic analysis on a data set for the 
diurnal and semidiurnal cycles as explained earlier (see equation (20)). Suff- 
icient data for 9 different data sets were found from MEN data. The results 
of the harmonic analysis compared very well with previous work that had been 
performed on essentially the same data^^ . 

After the tidal contributions had been extracted from a data set, the 
remaining data gave a series of irregular winds profiles as a function of time. 
The time intervals between the irregular data profiles provides the time separ- 
ation needed for the calculation of the time structure function. 

Analysis . The calculated time structure functions for the northward and 
eastward velocity components are shown in Figure 20. Both functions show a 
definite oscillatory motion which indicates the existence of significant motion 
with dominant periods. The dominant time scale (half period) is about 2-3 
hours. However, the oscillatory motion is centered on a value which is 
sufficiently high to indicate that there exists also irregular motion with a 
continuous distribution of periods. The lines drawn in the figure have slopes 
of 1.07 for the northward winds and 1.63 for the eastward winds. However these 
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values are highly uncertain due to the lack of data at short time displacements. 

A time structure function shown in Figure 21 was constructed from the corre- 

(91) 

lation data of Greenhow and Neufeld . This was done by using equation (9) 

and taking the rms irregular velocity to be 25 m/s, the value reported by Green- 

how and Neufeld. These results have a power law behavior at small time dis- 
placements with a 2/3 power of At. This result appears to contradict the MRN 

results but is in agreement with the horizontal structure functions since an 

application of Taylor’s hypothesis would say that time structure functions and 
horizontal structure functions should behave similarly. The data of Figure 21 
indicates a dominant time scale of about 3 hours, slightly higher than the MRN 
data, but the lengthening of the time scale is to be expected at the higher 
altitudes . 

(44) 

Justus has calculated the time structure function of irregular winds 
for the altitude range near 100 km. His results are given in Figure 22. At 
this higher altitude there seems to be less contribution from a single strongly 
dominant period in the irregular motions and more of a continuous distribution 
of periods. The slope of the linear increase is 1.5 which tends to confirm 
the higher of the slopes found at the lower altitude. The time scale at the 
higher altitudes is significantly increased and is of the order of 10 hours. 
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Figure 22. Time structure of chemical release horizontal winds previously ana- 

(44) 

lyzed by Justus using a deviation- from- the -me an technique. 





7. THE EFFECTIVE FREQUENCY SPECTRUM. 
ENCOUNTERED BY A SPACECRAFT 


The analysis of the structure functions presented in the previous sections 
indicate that the irregular wind, pressure, density, and temperature fields in 
the 50 - 200 km region can be considered as a more-or-less continuous spectrum, 
except that in certain cases, such as the vertical wave number spectrum in the 
45 - 65 km region, single modes of specific wavelength must be added to the 
otherwise continuous spectrum. This leads to the conclusion that the effects 
of the gravity wave field on spacecraft traversing these altitude regions can 
be analyzed in terms of the spectral techniques of turbulence. However, the 
results show that significant anisotropy exists in that the horizontal scales 
are significantly larger than vertical scales. In view of the nearly isotropic 
conditions in the horizontal, analysis in terms of axisymmetric spectra is 
suggested, with the axis of symmetry being in the vertical. 

The coordinates of the analysis are indicated in Figure 23. The x, y, z 
system has x and y axes in the horizontal plane and z in the vertical direction. 
The x direction has been conveniently chosen to be along the azimuth of the 
spacecraft trajectory. The elevation of the spacecraft trajectory is the 
angle 0 and the x f , y f , z 1 axes are those which are rotated through this angle 
0. Thus x 1 is along the direction of flight of the spacecraft. The displace- 
ment r along the x 1 axis is the distance of travel of the spacecraft in an 
arbitrary time interval x. With the speed of the spacecraft given by U, then 
r and t are related by r = Ux . The displacement r can be broken into com- 
ponent displacements r^ and r^ in the horizontal and vertical directions. 

Frequency Spectrum of Scalars . The frequency spectrum <f>(w) of a scalar 
property, say temperature T T , is related to the time correlation R(x) by the 
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Figure 23. Coordinate systems used in the analysis of observed frequency Spectra 
on spacecraft in a trajectory elevated at an angle 0. 


spectral transform relation 


0(u) 


<t ,2 > 


cos (JOT R(t) dx 


(35) 


The problem thus is to determine what the correlation R(x) is from such measur- 
able quantities as those presented in previous sections of this report. The 
assumption is made that the spacecraft speed U is sufficiently large that 
Taylor’s hypothesis can be used so that R(x) is the same as R(r) where r = Ux. 
It is further assumed that the correlation R(r) is the product of correlations 
which depend only on the components of r, i.e. 


R(r) = P 1 (^ 1 ) P 3 ( r 3 ) (36) 

If the horizontal and vertical structure functions D^(r^) and D^r^) have been 
determined instead of correlations, then, from equation (9), equation (36) can 
be written as 

R(r) = (1 - D 1 (r 1 )/2 <T' 2 >)(1 - D 3 (r 3 )/2<T' 2 >) (37) 

where r^ and r^ are given by 

r = U cos 0 x 

1 (38) 

r^ = U sin 0 x 

Thus any particular measured or model form for structure functions in (37) or 
correlation functions in (36) can be substituted into (35) . The effective fre- 
quency spectrum would, because of equation (38) , depend on the spacecraft 
velocity U and the elevation angle of the trajectory 0. It is beyond the 
scope of this report to evaluate actual frequency spectra. However, the re- 
sults of the previous sections indicate that a reasonable model would be given 
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by 

R(r) = [1 - (r 1 /r^ 0 ) m ] [1 - (r 3 /r 30 ) n ] (39) 

if < r^Q and r 3 < r 3 Q or R(r) = 0 otherwise, where r^ and r 3 Q are respect- 

* 

ively horizontal and vertical structure function scales and m and n are respect- 
ively the power law exponents of the horizontal and vertical structure functions 
at small scales. 

Frequency Spectrum of Velocity . Considerable additional complication is 
introduced when spectra of vector quantities, such as velocity, are considered. 

Not only must we consider longitudinal and transverse correlations but trans- 
formations between the horizontal - vertical system of coordinates (x, y, z) 
and the trajectory system of coordinates (x f , y ! , z f ). The spectrum of trans- 
verse (y ! ) velocity u^ is the same in either system of coordinates 

2 

» <u 2 > r 00 

<f>22 v — cos WT R 22^ t ^ dT (40) 


and the correlation function is likewise equal in both coordinates R 2 2 ? 
R 22 (t). But the "vertical" (z T ) velocity spectrum in the trajectory system 
must be transformed since we have 


♦ 33 ’(“> 




COS COT R^’Cx) 


dT 


(41) 


and as can be easily shown from the properties of coordinate rotation 


R33 ' ( T ) = C ° s2 6 R 33^ + sin2 6 R 11^ T ^ 

- sin 0 cos 0 [R^Ct) + R^ 3 (t)] 


(42) 


The correlations may also be expressed as function of r^ and r 3 given by equation 


66 



(38). However, the form of dependence required by the axisymmetric condition 

are somewhat more . complicated than (36). The form of the normalized correlation 

* (92) 

tensor (r) in an axisymmetric field is given by 

V? ’ Ar i r 3 + B 5 i3 s j3 + c s ij 


+ D( ri 6 j3 - rj 6 13 ) 


(43) 


where the axis of symmetry has been specified as the vertical, and the coeffic- 
ient functions A, B, C, and D are functions of the magnitude r and the vertical 
component magnitude |r 3 |. Of particular interest are the following specific 
terms of (43) 


R X1 (r) = A r 1 + C 

R 22 (r) = A r 2 + c 

R 33 (r) = A r ^ + B + C 

R l 3 (r) + R 31 (r) = 2A r 1 r 3 


(44) 


If we now assume that A = A^(r^) A 3 (r 3 ), B = B^(r^) B 3 (r 3 ) etc. then the 

easily measured correlations (and those which can be related by equation (9) 
to the structure functions discussed in the previous sections) can be written 
as 

2 . _ 


w 


+ R 22 (r 1 ) K 3 A l r 

w 

Rn(r 3 ) 

+ R 22 ( r 3 ) = 2c 3 

W 

R 33 (ri> 

= B 1 + C x 

V r 3> 

* E 33 (r 3 ) 

K l A 3 r 3 + B 3 + 


Where 


K = A x (0) and K 3 = A 3 (0) 


( 45 ) 
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Thus , B^, B^, C^, and are 6 unknowns required to determine the 

correlations and spectra completely, but the relations (45) , with the corre- 
* 

lations considered measured or otherwise known, constitute only 4 equations 
for these 6 unknowns. 

However, if we also assume that, in analogy with (36), the relations 

V r -> - KgCrp yr 3 )/2 

(46) 

(the denominator of 2 in the horizontal correlation is required to maintain 
the proper normalization ( 0) = 2 ) then it can be shown (see the Appendix) 

that 

Ai(r^) = Ry(r^) 

A 3 (r 3 ) = K 3 R H (r 3 )/2 

B 1 = B 3 = 0 (47) 

^l^ r l^ = ^V^ r l^ 

C 3 (r 3 } = V r 3 )/2 

where the constants and and microscales of the correlations (defined in 
the Appendix) must satisfy the relations 

K x K 3 = (2/X vl 2 - l/X H1 2 )/2 = 

(1/X h3 2 - 2/X v3 2 )/8 (48) 

Thus the correlation * s to employed in (40) for the deter- 

mination of the y T velocity spectrum is given by 
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r 22 (t) = . Ry(U cos 6 x) RjjCU Sin 0 t)/2 (49) 

and the correlations to be employed in (42) and (41) for the determination of the 

z * velocity spectrum are given by 

R n( x ) = (K-jK^ 2 + 1) RyC^) R R (r 3 )/2 

R 33 (t) = (K 1 K 3 r 3 2 + 1) R y (r 1 ) R R (r 3 )/2 (50) 

R 1 3 ( T ) + r 31 (t) = K^K 3 Ry(r 1 ) R H ( r 3 ) r 1 r 3 

where and are given by (38) . 

Thus any particular measured or model form of the correlation functions 
R^(r^), Ry(r^), and R^(r^) which are consistent in form with the equation 

(48) can be used in the above equations and substituted into (40) and (41) to 
determine the desired velocity spectra. The effective frequency spectra are, 
because of equation (38) , dependent on the spacecraft velocity U and elevation 
angle of trajectory 6. If horizontal and vertical structure functions have been 
determined instead of correlation functions, then from equation (9) the follow- 
ing relations can be employed instead: 

R H (r 1 ) = 1 - D R ( ri ) / [2 (<u 2 > + <v 2 >)] 

R R (r 3 ) = 1 - D R (r 3 ) /[2(<u 2 > + <v 2 >)] 

2 (51) 

R y (r 1 ) = 1 - D v (r x ) /2 <w > 

R y (r 3 ) = 1 - D y (r 3 ) /2 <w 2 > 

Again, it is beyond the scope of this report to evaluate actual frequency spectra. 
However* model correlations similar to equation (39) would be reasonable forms 
for the rough approximation of the measured results of the previous sections. 
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APPENDIX 


Under the assumption of equation (46) the normalized . correlations R^ and 


R v become, after substitution from (45) 


Rr(-) K 3 A 1 C 3 r l + 2C 1 C 3 


R v (r) = K x A 3 (B 1 + C 1 ) r 3 + B 1 B 3 + ^ C 3 


(Al) 


+ B 1 C 3 + Cj B 3 


In order for these relations to be equivalent to the defining relations for R^ 
and Ry 


R R (r) = A r R 2 + 2C 
S v (r) = A r 3 2 + B + C 


(A2) 


where, because of the coordinates employed in Figure 23, = r^, then because 

we have assumed A * A^, B = B-^ B^, anc ^ C = C^, the following relations 

apply between the coefficient function 


A 3 * K 3 C 3 


Aj - + 03) 


B 1 C 3 + C 1 B 3 - 0 


(A3) 


Since was defined as A^(0) anc ^ ^3 as A ^(0) then (A3) shows that C^(0) = 1 
and (0) + 0^(0) = The normalization R^(0) = 2 together with the value 

C^(0) = 1 implies that C^(0) = 1. Hence B^(0) = 0. The last relation of 

(A3) may be written as 
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Bprp B 3 ( r 3 ) 

C l^ r l^ C 3 ^ r 3> 


constant 


(A4) 


where the constancy is required by the fact that the first term depends only on 
r^ and the second depends only on r^. Now since (A4) must hold for all values 
of r^ then the value of the constant is B^(0) / C^(Q) which is zero. Thus 


B 1 (r 1 ) = 0 

B 3 (r 3 ) = 0 


(A5) 


for all values of and r^, and the first two relations of (A3) thus become 


A l (r l> = K x C^rp 


A 3 ( r 3 ) = K 3 C 3 (r 3 ) 


(A6) 


and (45) simplifies to 


V r i } 

= (Kj_ K 3 r x + 2) ^(rp 

(A7) 

W 

= 2 C 3 (r 3 ) 

(A8) 

V r i> 

= Cprp 

(A9) 


- (K x K 3 r 3 2 + 1) C 3 (r 3 ) 

(A1Q) 


Relations (A5) , (A6) , (A8) , and (A9) thus give the coefficient functions A^, A 3 , 
B^, B 3> C^, and C 3 in terms of the correlations R^(r 3 ) and R^(rp . Second order 
differentiation of the ratio of (A7) to (A9) and the ratio of (A10) to (A8) 
establish two equations which relate the constants and K 3 to the derivatives 
of the correlation functions. Since the resulting equations must be valid at 
all values of r^ and r 3 then they can be evaluated by setting r^ and r 3 to zero 


71 



and by noting that R^( 0) - 1, R^CO) = 2, and all first order derivatives of 

at r^ = 0 or r^ = 0 are zero. The result- 
ing equations are given by (42) where the second derivative terms evaluated at 
zero are defined in terms of microscales. 


and with respect to or r^ 


(d 2 R H /dr 1 2 ) 

(d 2 \/dr 3 Z ) 
(d 2 R v /dr 1 2 ) 
(d 2 R v /dr 3 2 ) 







(All) 
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